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F. Yusef-Zadeh 1 , M. Wardle, 2 , H. Bushouse 3 , CD. Dowell 4 , and D. A. Roberts 1 

ABSTRACT 

We have investigated the nature of flare emission from Sgr A* during multi- 
wavelength observations of this source that took place in 2004, 2005 and 2006. We 
present evidence for dimming of submm and radio flux during the peak of near-IR 
flares. This suggests that the variability of Sgr A* across its wavelength spectrum is 
phenomenologically related. The model explaining this new behavior of flare activity 
could be consistent with adiabatically cooling plasma blobs that are expanding but 
also partially eclipsing the background quiescent emission from Sgr A*. When a 
flare is launched, the plasma blob is most compact and is brightest in the optically 
thin regime whereas the emission in radio/submm wavelengths has a higher opacity. 
Absorption in the observed light curve of Sgr A* at radio/submm flux is due to the 
combined effects of lower brightness temperature of plasma blobs with respect to 
the quiescent brightness temperature and high opacity of plasma blobs. This implies 
that plasma blobs are mainly placed in the magnetosphere of a disk-like flow or 
further out in the flow. The depth of the absorption being larger in submm than in 
radio wavelengths implies that the intrinsic size of the quiescent emission increases 
with increasing wavelength which is consistent with previous size measurements of 
Sgr A*. Lastly, we believe that occultation of the quiescent emission of Sgr A* at 
radio/submm by IR flares can be used as a powerful tool to identify flare activity at 
its earliest phase of its evolution. 

Subject headings: accretion, accretion disks — black hole physics — Galaxy: center 



1. Introduction 

Observations of stellar orbits in the proximity of the enigmatic radio source Sgr A* located at 
the very dynamical center of our galaxy have established that it is a 4 xlO 6 M black hole 
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sists of quasi-steady and variable components (e.g. 



2009 ). The emission from Sgr A* con 
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Yusef-Zadeh et al.ll2010l ). The variable component is now detected in almost all wavelength bands 



and is used to estimate the physical quantities of the gas flow, the radiation mechanism in different 
wavelength bands and the scale lengths at which variable emission operates. One of the puzzling 
aspect of the flaring activity of Sgr A* is the physical nature of the variable emission in different 
wavelengths and how its spectrum evolves with time. Furthermore, it is not clear whether flares 
are deeply embedded in the accretion disk of Sgr A* or are physica l ly dis t inct from the region 



where the bulk of the emission aris es (|Yusef-Zadeh et al.l l2006al . 120081 . |2009| ; iMarrone et al.l 120081 ; 



Eckart et al.ll2008l : iTrap et al.ll2010h . 



Flare emission at near-IR (IR) wavelengths is due to optically thin synchrotron emission and the 
submm frequency is considered to be the dividing line between optically thin emission at high 
frequencies (near-IR), and optically thick emission at lower frequencies (submm/radio). To examine 
the possibility that there is a causal association between radio/submm variability and IR flaring 
activity, we studied in detail the light curves of flare emission from optically thick and thin plasma 
where there was sufficient overlap between IR/X-ray and submm/radio wavelengths. We searched 
for signatures of optical depth effect against the background quiescent emission at the earliest phase 
of the evolution of flares. 



Data Reduction and Results 



The primary purpose of our archival search was to find simultaneous observations of Sgr A* at 
submm, radio and IR wavelengths. We found three experiments that had overlaps between these 
wavelengths bands. Almost all the data presented here have already been published previously. 

One multi- wavelength observing campaign was carried out on 2006 July 17 during which there was 
wavelength coverage between the Keck, Chandra, the Caltech Submillimeter Observatory (CSO), 
the Submillimeter Array (SMA), and the Very Large Array (VLA) o f the National Radio Astr onomy 
Observatorji (NRAO) (jHornstein et alJl2007l : iMarrone et allhoosl : lYusef-Zadeh etaD bood ) . The 
CSO data are reanalyzed with a higher time sampling and the VLA data are recalibrated. The 
light curves presented at 7 and 13mm have selected data corresponding to uv spacings greater than 
lOOkA. 

The second multi-wavelength observations took place on 2005 July 31 employing the CSO, the 
SMA, the Keck II and Very Large Telescopes (VLT). The time coverage on this day was excellent 
as there was a total of ten continuous hours of observations wi t h the Keck I I and VLT telescopes 
at IR wavelengths (jHornstein et al.l 120071 : IMarrone et al.1 120081 : iMeyer et al.l 120081 ) . We searched 
for archival radio data from the VLA and found 7mm (43GHz) and 13mm (12GHz) data taken 
for a temporal coverage of about four hours on 2005, July 31. These data were taken in the C 
configuration using fast switching technique and a phase calibrators (0173+331, 17444-31166) to 
take out rapid atmospheric fluctuations. The light curves presented here are self-calibrated in phase 
using data corresponding to uv spacings > 50kA. 



Lastly, the campaign that took place on 2004 September 4 made simultaneous observations with 
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the Hubble Space Tele scope (HST/NICMOS) an d CSO. The submm data are reanalyzed with a 
better time sampling (jYusef-Zadeh et al.l l200 6b). The NICMOS data include photometry from 
three filters - F160W, F187N, and F190N - which are centered at wavelengths of 1.6, 1.87, and 
1.90/im. The aperture photometry measurements have been corrected to an absolute flux scale, 
including an extinction correction that assumes A^-=3.3 magnitudes for Sgr A*. 

Three examples in which there is evidence for anti-correlation of optically thin and thick flare 
emission are presented. The first set of light cu rves show radio observations of flaring activity that 
took place on July 31, 2005. iMeyer et al.l (J2008|) presented a 600min light curve by combining Keck 
II and VLT data. Sgr A* was highly active during this period showing several IR flares. Figure la 
shows IR and radio light curves. The bottom panels show individual 7mm and 13mm light curves 
in two different colors. Radio measurements indicate a remarkable trend that the emission is higher 
when there is a low-level flaring activity and is lower when there is a strong IR flare. This trend 
is best noted at times less than 5h UT and greater than 6h UT. The flux changes by about 150 
and 100 mJy at 7 and 13mm wavelengths, respectively. A more detailed view of flaring activity 
between 5.5h and 7.5h UT is shown in Figure lb. The light curves cit cit 2.2 ^m, 7 and 13mm are 
presented with sampling time of 60s, 30sec and 30sec, respectively. The anti-correlation is evident 
clearly, first by the rise and fall of IR and radio activity during 5.5 and 7.5h UT. A more detailed 
examination of the radio emission anti-correlating with two of the strongest IR flares can best be 
viewed near 6h UT. Radio flux rises between the two peaks of flare emission at 5.8h and 6.2 UT. 
The fast fluctuations in the flux of Sgr A* at radio wavelengths appears to track the variation of 
flux at IR wavelengths. To clarify this relationship better, radio flux of Sgr A* are reversed before 
the change in radio flux is presented in Figure lc. The "correlation" of radio and IR flux is shown 
in remarkable clarity, given that the quiescent flux and the net variable flux are different at 7 and 
13mm wavelengths. The apparent discrepancy between 7mm and 2.2/iin variation near 7.1h UT 
may be due to poor time sampling ~ 3min of Keck data and a sudden change of observations 
from VLT to Keck which starts around 6.8h UT. The anti-correlation between near-IR and 1.3mm 
data (Marrone et al. 2008) was not as clear as that of radio and near-IR data but 1.3mm data 
could suffer from poor sampling and complication from time delayed emission. To examine the 
anti-correlation statistically, Figure 2a, b show the cross correlation of the IR data with 7mm and 
the merge d 7 and 13mm da ta, respectively. The cross-correlation analysis use the ZDC function 
algorithm ([Alexander! 119971 ). A minimum in the likelihood value is identified at a zero time lag, 
strengthening the anti-correlation of IR and radio data. 

Th e second evidence for dimming of optically thick light comes from observations made in 2006 July 
17 (jMarrone et al.l 120081 ; lYusef-Zadeh et al.ll2008l ). The authors showed evidence of a time delay 
between the peaks of IR/X-ray emission with respect to radio/submm emission. Figure 3a shows 



simu ltaneous light curves of Sgr A* at X-rays between 2 and 8 keV, IR K' band (jHornstein et al 



20071 ). 850/iin, 1.3mm, and 7mm. In these observations, as noted in the top panel of Figure 3a, 



the second half of the X-ray flare was observed with Keck and showed evidence of flaring at IR 
wavelengths. This figure shows a dip in the flux of mm and submm emission, during which an X-ray 
flare is detected. A close-up view during the peak emission at X-rays between 5:30 and 8h UT is 
shown in Figure 3b. The light curves at 1.3mm and 850/im track each other well. However, we 
note the 7mm light curve seems to be somewhat shifted in time by with respect to the 1.3mm light 
curve. The poor sampling of 1.3mm data between 6 and 7h UT makes the comparison difficult. The 
largest drop in the 1.3mm flux is noted at 6.7h UT. The average dimming at 850/im, 1.3mm and 
7mm range between 0.5-1 Jy, 0.5 Jy and 200-300 mJy, respectively. The 13mm data also showed a 
very similar light curve to that of 7mm with a flux variation of about 100-150 mJy. 
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Although we note a general trend in the anti-correlation of X-ray and optically thick emission, this 
anti-correlation should have been evident between IR and optically thick emission. This is because 
X-ray emission is processed IR emission and does not exactly tract the flux of a IR flare. The 



dura t ion of a typical X-ra y flare is generally shorter than in IR wavelengths (jDodds-Eden et al 



2009; iMarrone et al.ll2008l ). The 2.2//m flare of Figure 3 indicates the duration of flare is longer 



than that observed in X-rays by at least 30 minutes on the decaying side of the flare. Unfortunately, 
there is no IR data available at times earlier than 6:45h UT. Nevertheless, the rise of flux at 7mm, 
1.3mm and 850/um between 6h and 6:45h UT is consistent with the decline in the IR flux. Using 
the anti-correlation between IR and radio flux, we expect that the IR flare should have started 
around 5h UT with a flare duration of about two hours. Another feature we note in Figure 3 is a 
possible trend that the dimming of light appears to end first at longer wavelengths. However, this 
needs to be confirmed as SMA data are not sampled well during the dip in the quiescent flux of 
Sgr A*. 

Lastly, the light curves showing evidence of dimming at submm wavelengths is also detected in an 
observing campaign that took place on 2004, September 4. The light curves presented in Figure 
4a show emission at 850/um and three IR wavelengths at 1.6, 1.87 and 1.90/im. Two strong IR 
flares and on e submm flare are d etected near 4.2h UT, 7.6h UT and 7h UT, respectively, (see also 



Figure 10 of lYusef-Zadeh et al.l (l2006bl ). A previous analysis of submm data using long sampling 
time of 20min in the construction of the light curve suggested the possibility that the submm flare 
is the counterpart either to flare at 4.3h UT delayed by ^200min with respect to the first IR flare 
or to the second IR flare at 7.6 with no time delay (jYusef-Zadeh et al.l l2006bl ) . However, the new 
analysis shows a new feature that had not been recognized before. The 850/iin light curve indicates 
a decrease in the quiescent flux of Sgr A* by ~ 500 mJy during which the second IR flare occurred. 
We also note a decrease in the submm flux at the start of IR observations when there is a detection 
of a flare. Figure 4b shows detailed views of the light curves to demonstrate the brightening and 
dimming of IR and submm flux during flaring activity near 7.6h UT, respectively. This figure 
provides another support for a general trend in the anti-correlation between IR flare emission and 
submm flux of Sgr A*. Exceptions to this trend are two data points at 7.8h UT during which 
submm flux is lower than the peak IR flux near 7.5h UT. Unfortunately, the data are not sampled 
well at these wavelengths. Another effect that may complicate the analysis of absorption in light 
curves at radio and submm wavelengths could be due to flaring events that are time delayed and 
are superimposed on an absorbing feature. 



3. Discussion 

The infrared wavelength band is the most important part of the spectrum of flare emission from 
Sgr A*. This is because the flare is first launched in IR wavelengths followed by radio and submm 
emission being delayed with respect to the peak in IR emission. In this model, flaring at a 
given frequency is produced through the adiabatic expansion of an initially optically thick blob 
of synchrotron-emitting relativistic electrons. The expanding hot plasma model of Sgr A* flares 
has explained the nature of time delay of flare emission in the optically thick regime. However, 
the relationship between IR and radio/submm flare emission has not been well established. Recent 
work has suggested that that the optically thin component of flare emission is unrelated to the 
evolution of optically thick e mission and that t hat IR and radio/submm flares emission operate 



independent of each other (e.g. lSabha et al.ll2009l ). The light curves presented here demonstrate the 
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dimming of the quiescent flux of Sgr A* at radio/submm wavelengths while IR flares are observed. 
The new measurements show that the observed flares in different frequencies operate under the 
same physical mechanism. 

A likely scenario for the drop in flux seen at radio and submm frequencies during IR flaring is that 
the electron population responsible for the flaring is optically thick at lower frequencies, partially 
obscures the submm-radio emission from the accretion flow, and has a lower brightness temperature. 
To explore this picture we construct a simple two-component model. The first component is a 
homogeneous, foreground, synchrotron emission region responsible for the IR flare. We characterize 
this by an E~ 2 electron energy spectrum of relativistic electrons, a magnetic field strength B, and 
radius R. This component is responsible for the IR flare and also for producing the dips in the 
emission from Sgr A*. The second, background, component is responsible for the quiescent flux 
from Sgr A* in the radio and submm, is also homogeneous and is parametrized by its flux density 
S qu and radius R qu , both of which are frequency-dependent. For simplicity, we assume that the 
two components are concentric, so that the fraction of the quiescent background source covered by 
the flaring component is / = (R/R qu ) 2 when R < R qv and / = 1 otherwise. 

The net flux received from Sgr A* in this simple model consists of contributions from the flaring 
region, and the obscured and unobscured fractions / and 1 — / of the quiescent source: 

S u = £IJ U (I - exp(-r„)) + (1 - f)S qv + fSq V exp(-T„) (1) 

= S qu + (tlJ„ - fS qv ) (1 - exp(-7v)) , (2) 

where the flare region has source function J u , optical depth t u , and subtends solid angle Q = ttR? jd 2 
at Earth, where d = 8 kpc is our adopted distance to the Galactic Center. 

The synchrotron source function for an electron number density spectrum n(E) oc E~ 2 is 

P- 

where 

_ / A-Km e cu \ 1//2 o , A s 

E " = hsr) m - c (4) 

is the typical energy of the electrons responsib le for synchrotron emission at frequency is, and the 



J u = 0.646 % E u , (3) 

cr 



other symbols have their usual meaning (e.g. iRybicki fc Lightmanl Il986l ). We write the optical 



depth as t u = k u R, where the synchrotron absorption coefficient is 

2 e 3 B 



9 m e u 2 



n{E v ) . (5) 



Note that we have implicitly set the electron pitch angle to 90 degrees in equations (jSJ) and 
Finally, we adopt a fixed ratio between the energy density of the electron population, i.e. U e = 
f E 2 En(E) dE, and the energy density of the magnetic field, i.e. Ub = B 2 /8tt, so that 

between energies E\ and E2, where 4> = (Us/Ue) ln(£ , 2/ 'Ex). 

The flux in the infrared is dominated by optically thin flare emission, and equation ([2]) reduces to 

S u = QJ u t u , (7) 
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This allows us to use equations (J3j) — dSJ) to compute the flare region's magnetic field strength, optical 
depth and source-function flux 

B = 38(^ 3 5 10 ) 2/7 R u 6/7 G, (8) 
t v = 0.lU(cf> 3 )-V 7 sV 7 Ru 17/7 v^, (9) 

nj u = i.7(^s 10 )- 1/7 Rll /7 4io Jy. (io) 

where we have scaled the results using R12 = -R/10 12 cm, ^3 = 0/3, S\o = Sv/10mJy at 2.2 //m, 
and ^350 = ^/350GHz. The weak dependence on the equipartition factor <j) mean that our results 
are insensitive to the ratio of the energy densities in the magnetic field and electron population. 

Having determined the properties of the flare source region, we can compute its obscuring effect 
on the radio-submm quiescent emission using (j2J). By way of illustration we adopt = 3, lOmJy 
flux of at 2.2 fj,m, and quiescent fluxes S qu = 3.0 and 1.8 Jy at v = 350 and 43 GHz, respectively. 
At each frequency, for a given choice of quiescent source radius R qu we compute the flare region 
radius R that would produce a specified drop in the total flux. The results are plotted in Fig. 5 for 
flare-induced drops of 0.05, 0.1, and 0.2 Jy at 350 GHz (850//m) and 43 GHz (7mm). 

At 350 GHz, absorption of ~ 0.1 Jy during an IR flare implies that both R and R qv ~ 10 12 cm. 
The latter is consistent with the size of Sgr A* inferred from VL BI measurements at 230 GHz, 



which found a major axis FWHM ~ 37 //as (jDoeleman et al.ll2008l ) and in conjunction with lower 



frequency measurements ( Bower et al.ll2006l ) infer a z/ 1 - 44 ^ 007 wavelength dependence, implying a 
major axis FWHM of ~ 20 //as at 350 GHz, equivalent to 2.4 x 10 12 cm. While choices of smaller 
R in principle produce similar absorption levels, the magnetic field strength scales as it! -6 / 7 , and 
exceeds 100 G for R < 3 x 10 11 cm. 

The 43 GHz curves in Figure 5 imply that values of R much below 10 12 cm require unrealistically 
small size for Sgr A* at 43 GHz. The radius of a flare R ~ 10 12 cm , on the other hand, implies that 
R qu « 5 x 10 12 cm at 43 GHz. While this is several times smaller than the measured semi major 
axis, ~ 2 x 10 13 cm, the semi-minor axis is not yet determined and may be a few times less, so this 
discrepancy is not as large as may first appear. 

This simple model shows that occultation of Sgr A* by an IR flare source can produce noticeable 
dips in the flux from Sgr A* at radio and sub- mm frequencies. The major uncertainties are the 
electron spectral index, which is used to infer the optical depth at radio and submm from the IR 
flux density, and the requirement that the IR flare region must overlie the radio emission. Note 
however that for a given IR flux density, a steeper electron spectrum implies greater optical depth 
at low frequencies. This latter would more naturally occur if the flare emission occurred in a 
magnetosphere of a disk-like flow or further out in the flow. One would expect that the geometry 
would not be favorable all of the time. 



A rec ent model has considered time dependent model of flare emission from Sgr A* (IDodds-Eden et al 



2010) . They argue that flares are generated from the d issipation of the m agnetic field due to the 



reconnection of the field in the accretion disk (see also lYuan et al.1 12009). This model predicts a 
dimming of flux at submm flux. Although their picture is quite consistent with some of the light 
curves presented here, more detailed observations are required to distinguish between these two 
models. In other words, the question is whether the occultation is to due to an absorbing blob of 
hot plasma against a background quiescent flux of Sgr A* or due to a decrease in the flux as a 
a consequence of the dissipation of the magnetic field due to reconnection of magnetic field lines. 
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Future detailed modeling and observations with long time coverage should be able to address the 
difference between these two models. 
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Fig. 1. — : (a) Top The IR data was extracted from Figure 3 of lMever et al.l (|2008l ). The sampling 
time of 7 and 13mm data is 300sec. (b) Bottom Left The same as in (a) except that a close-up 
view is presented. The sampling time for the IR, 7mm and 13mm are about lmin, 30sec and 
30sec, respectively, (c) Bottom Right Similar to (b) except that radio data are reversed by first 
multiplying by -1 and then adding a flux of 1.8 and 1.3 at 7 and 13mm, respectively. 
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Fig. 2. — : (a) Left The cross correlation of IR and 7mm light curves is shown in the bottom panel. 
(b) Right Similar to (a) except that radio data has merged the 7 and 13mm data by adding a 
constant flux of 0.4 Jy to the 13mm data. 
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Fig. 3. — : (a) Left The light curves of Sgr A* is shown in five different wavelengths The X-ray (2-8 
keV), IR, submm (850 /i m), 1.3mm and 7 mm data are binned every 207sec, 148sec, 120sec, 2 50sec 
and 250sec, respectively (jHornstein et all 120071 ; Marrone et al-lbood ; lYusef-Zadeh et aliboosl ). (b) 
Right The same as (a) except that a close up view of the light curves including a IR data at 1/ 
band. 




UT Time (Hours) UT Time (Hours) 

Fig. 4.—: (a) Left The simultaneous IR and 850/um light curves using NICMOS/HST and CSO 
on 2004, September 4. The 1.6, 1.87 and 1.90, 850//m sampling time is 32sec, 192sec, 192sec and 
280sec, respectively, (b) Right Similar to (a) except that a close-up view is shown. 




Fig. 5. — : Combination of IR flare region radius and Sgr A* quiescent emission radius needed to 
produce drops in flux density at 43 and 350 GHz of 50, 100 and 200 mJy (short-dashed, solid, and 
long-dashed curves, respectively.) during an infrared flare. 



